On the Use of Correlation and MI as a Measure of Metabolite—Metabolite Association for Network Differential Connectivity Analysis

Metabolite differential connectivity analysis has been successful in investigating potentialmolecular mechanisms underlying different conditions in biological systems. Correlation and MutualInformation (MI) are two of the most common measures to quantify association and for buildingmetabolite—metabolite association networks and to calculate differential connectivity. In this study,we investigated the performance of correlation and MI to identify significantly differentially connected metabolites. These association measures were compared on (i) 23 publicly available metabolomic data sets and 7 data sets from other fields, (ii) simulated data with known correlation structures,and (iii) data generated using a dynamic metabolic model to simulate real-life observed metabolite concentration profiles. In all cases, we found more differentially connected metabolites when using correlation indices as a measure for association than MI.We also observed that different MI estimation algorithms resulted in difference in performance when applied to data generated using a dynamic model. We concluded that there is no significant benefit in using MI as a replacement for standard Pearson’s or Spearman’s correlation when the application is to quantify and detect differentially connected metabolites

Approaches to Sample Size Determination for Multivariate Data:Applications to PCA and PLS-DA of Omics Data

Sample size determination is a fundamental step in the design of experiments. Methods for sample size determination are abundant for univariate analysis methods, but scarce in the multivariate case. Omics data are multivariate in nature and are commonly investigated using multivariate statistical methods, such as principal component analysis (PCA) and partial least-squares discriminant analysis (PLS-DA). No simple approaches to sample size determination exist for PCA and PLS-DA. In this paper we will introduce important concepts and offer strategies for (minimally) required sample size estimation when planning experiments to be analyzed using PCA and/or PLS-DA.</p

Considering Horn’s parallel analysis from a random matrix theory point of view

Horn’s parallel analysis is a widely used method for assessing the number of principal components and common factors. We discuss the theoretical foundations of parallel analysis for principal components based on a covariance matrix by making use of arguments from random matrix theory. In particular, we show that (i) for the first component, parallel analysis is an inferential method equivalent to the Tracy–Widom test, (ii) its use to test high-order eigenvalues is equivalent to the use of the joint distribution of the eigenvalues, and thus should be discouraged, and (iii) a formal test for higher-order components can be obtained based on a Tracy–Widom approximation. We illustrate the performance of the two testing procedures using simulated data generated under both a principal component model and a common factors model. For the principal component model, the Tracy–Widom test performs consistently in all conditions, while parallel analysis shows unpredictable behavior for higher-order components. For the common factor model, including major and minor factors, both procedures are heuristic approaches, with variable performance. We conclude that the Tracy–Widom procedure is preferred over parallel analysis for statistically testing the number of principal components based on a covariance matrix.<p>Horn’s parallel analysis is a widely used method for assessing the number of principal components and common factors. We discuss the theoretical foundations of parallel analysis for principal components based on a covariance matrix by making use of arguments from random matrix theory. In particular, we show that (i) for the first component, parallel analysis is an inferential method equivalent to the Tracy–Widom test, (ii) its use to test high-order eigenvalues is equivalent to the use of the joint distribution of the eigenvalues, and thus should be discouraged, and (iii) a formal test for higher-order components can be obtained based on a Tracy–Widom approximation. We illustrate the performance of the two testing procedures using simulated data generated under both a principal component model and a common factors model. For the principal component model, the Tracy–Widom test performs consistently in all conditions, while parallel analysis shows unpredictable behavior for higher-order components. For the common factor model, including major and minor factors, both procedures are heuristic approaches, with variable performance. We conclude that the Tracy–Widom procedure is preferred over parallel analysis for statistically testing the number of principal components based on a covariance matrix.</p

A Metabolomic Perspective on Coeliac Disease

Metabolomics is an “omic” science that is now emerging with the purpose of elaborating a comprehensive analysis of the metabolome, which is the complete set of metabolites (i.e., small molecules intermediates) in an organism, tissue, cell, or biofluid. In the past decade, metabolomics has already proved to be useful for the characterization of several pathological conditions and offers promises as a clinical tool. A metabolomics investigation of coeliac disease (CD) revealed that a metabolic fingerprint for CD can be defined, which accounts for three different but complementary components: malabsorption, energy metabolism, and alterations in gut microflora and/or intestinal permeability. In this review, we will discuss the major advancements in metabolomics of CD, in particular with respect to the role of gut microbiome and energy metabolis

SAPP: functional genome annotation and analysis through a semantic framework using FAIR principles

There are currently more than 150.000 sequenced genomes available from which considerable amounts of information can be extracted. However, annotation information is often not interoperable, static, lacks provenance and is quickly outdated. Keeping these datasets up-to-date, and interoperable is a challenging ..

Strategies for individual phenotyping of linoleic and arachidonic Acid metabolism using an oral glucose tolerance test

The ability to restore homeostasis upon environmental challenges has been proposed as a measure for health. Metabolic profiling of plasma samples during the challenge response phase should offer a profound view on the flexibility of a phenotype to cope with daily stressors. Current data modeling approaches, however, struggle to extract biological descriptors from time-resolved metabolite profiles that are able to discriminate between different phenotypes. Thus, for the case of oxylipin responses in plasma upon an oral glucose tolerance test we developed a modeling approach that incorporates a priori biological pathway knowledge. The degradation pathways of arachidonic and linoleic acids were modeled using a regression model based on a pseudo-steady-state approximated model, resulting in a parameter A that summarizes the relative enzymatic activity in these pathways. Analysis of the phenotypic parameters As suggests that different phenotypes can be discriminated according to preferred relative activity of the arachidonic and linoleic pathway. Correlation analysis shows that there is little or no competition between the arachidonic and linoleic acid pathways, although they share the same enzyme

Group‑wise ANOVA simultaneous component analysis for designed omics experiments

Modern omics experiments pertain not only to the measurement of many variables but also follow complex experimental designs where many factors are manipulated at the same time. This data can be conveniently analyzed using multivariate tools like ANOVA-simultaneous component analysis (ASCA) which allows interpretation of the variation induced by the different factors in a principal component analysis fashion. However, while in general only a subset of the measured variables may be related to the problem studied, all variables contribute to the final model and this may hamper interpretatio

Necrotizing soft tissue infections - a multicentre, prospective observational study (INFECT) : Protocol and statistical analysis plan

Background: The INFECT project aims to advance our understanding of the pathophysiological mechanisms in necrotizing soft tissue infections (NSTIs). The INFECT observational study is part of the INFECT project with the aim of studying the clinical profile of patients with NSTIs and correlating these to patient-important outcomes. With this protocol and statistical analysis plan we describe the methods used to obtain data and the details of the planned analyses. Methods: The INFECT study is a multicentre, prospective observational cohort study. Patients with NSTIs are enrolled in five Scandinavian hospitals, which are all referral centres for NSTIs. The primary outcomes are the descriptive variables of the patients. Secondary outcomes include identification of factors associated with 90-day mortality and amputation; associations between affected body part, maximum skin defect and Laboratory Risk Indicator for Necrotizing Fasciitis (LRINEC) score and 90-day mortality; 90-day mortality in patients with and without acute kidney injury (AKI) and LRINEC score of six and above or below six; and association between affected body part at arrival and microbiological findings. Exploratory outcomes include univariate analyses of baseline characteristics associations with 90-day mortality. The statistical analyses will be conducted in accordance with the predefined statistical analysis plan. Conclusion: Necrotizing soft tissue infections result in severe morbidity and mortality. The INFECT study will be the largest prospective study in patients with NSTIs to date and will provide important data for clinicians, researchers and policy makers on the characteristics and outcomes of these patients.</p

Of Monkeys and Men: A Metabolomic Analysis of Static and Dynamic Urinary Metabolic Phenotypes in Two Species

Background Metabolomics has attracted the interest of the medical community for its potential in predicting early derangements from a healthy to a diseased metabolic phenotype. One key issue is the diversity observed in metabolic profiles of different healthy individuals, commonly attributed to the variation of intrinsic (such as (epi)genetic variation, gut microbiota, etc.) and extrinsic factors (such as dietary habits, life-style and environmental conditions). Understanding the relative contributions of these factors is essential to establish the robustness of the healthy individual metabolic phenotype. Methods To assess the relative contribution of intrinsic and extrinsic factors we compared multilevel analysis results obtained from subjects of Homo sapiens and Macaca mulatta, the latter kept in a controlled environment with a standardized diet by making use of previously published data and results. Results We observed similarities for the two species and found the diversity of urinary metabolic phenotypes as identified by nuclear magnetic resonance (NMR) spectroscopy could be ascribed to the complex interplay of intrinsic factors and, to a lesser extent, of extrinsic factors in particular minimizing the role played by diet in shaping the metabolic phenotype. Moreover, we show that despite the standardization of diet as the most relevant extrinsic factor, a clear individual and discriminative metabolic fingerprint also exists for monkeys. We investigate the metabolic phenotype both at the static (i.e., at the level of the average metabolite concentration) and at the dynamic level (i.e., concerning their variation over time), and we show that these two components sum up to the overall phenotype with different relative contributions of about 1/4 and 3/4, respectively, for both species. Finally, we show that the great degree diversity observed in the urinary metabolic phenotype of both species can be attributed to differences in both the static and dynamic part of their phenotyp